An irreversible magnetic-field dependence of low-temperature heat transport of 
spin-ice compound Dy 2 Ti 2 7 in a [111] field 
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We study the low-temperature thermal conductivity (k) of Dy2Ti2C>7 along and perpendicular 
to the (111) plane and under the magnetic field along the [111] direction. Besides the step-like 
decreases of k, at the field-induced transitions from the spin-ice state to the kagome-ice state and 
then to the polarized state, an abnormal phenomenon is that the k,(H) isotherms show a clear irre- 
versibility at very low temperatures upon sweeping magnetic field up and down. This phenomenon 
surprisingly has no correspondence with the well-known magnetization hysteresis. Possible origins 
for this irreversibility are discussed; in particular, a pinning effect of magnetic monopoles in spin ice 
compound by the weak disorders is proposed. 



PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 

I. INTRODUCTION 

The isolated magnetic charges or the magnetic 
monopoles, in contrast to their electrical counterparts, 
have not been observed in nature. A recent discovery 
was that the magnetic monopoles could emerge from the 
collective excitations of spin ice* 1 - - — Although this kind 
of magnetic-monopole excitations is not completely the 
same as the free elementary particles, this finding is in- 
teresting and very important in the sense that it was the 
first time to signify the magnetic monopoles in real space. 

E>y2Ti207 (DTO), as a spin ice material, has already 
become the focus of condensed matter physics 1 - and re- 
ceived a lot of attentions during past ten years due to its 
fantastic physical properties.— ~— DTO belongs to the 
family of pyrochlore titanates with the Dy 3+ spins locat- 
ing at the vertexes of tetrahedra, which consist of trian- 
gular and kagome planes stacked alternatingly along the 
[111] direction. At low temperatures, the Dy 3+ spins are 
Ising anisotropic and form a macroscopically degenerated 
"2-in, 2-out" spin-ice ground state, which is attributed to 
the frustration effect between the nearest-neighboring an- 
tiferromagnetic interaction and the ferromagnetic dipolar 
interaction^ - — Once the flipping of a spin occurs, a local 
"3-in, 1-out" or "3-out, 1-in" spin configuration forms, 
which is equivalent to yielding two opposite magnetic 
monopoles in the adjacent tetrahedraj^r— as shown in 
Fig. 1(a). The monopole pairs can be separated and dif- 
fuse freely in the spin-ice state, which can be viewed as 
a "vacuum" free of charges by continuous spin flip and 
behaves a random walk process in zero field* 1 ^ A partic- 
ular case of interest is with the magnetic field along the 
three- fold axis or the [111] direction. With increasing 
field, there are two successive transitions from the spin- 
ice state to the kagome-ice state and then to the fully- 
polarized stated - — In the kagome-ice state, the mag- 




FIG. 1: (Color online) (a) "2-in, 2-out" spin configuration of 
the spin-ice ground state in Dy2Ti207. The spin-easy axis is 
along the three- fold axis or the [111] direction. Two neigh- 
boring tetrahedra are shown and the arrows stand for Dy 3+ 
spins on the vertexes. When one spin is flipped, indicated by 
the open arrow, by thermal fluctuations or external magnetic 
field, a pair of magnetic monopoles are formed and shown by 
the red and blue dots, (b) In magnetic field perpendicular 
to the (111) plane, the magnetic monopoles can move in the 
plane by successive spin flipping. The dashed red and blue 
arrows, accompanied with the dashed open circles, show the 
routines of positive and negative monopoles. 

netic field pins the spins on the triangular planes along 
the field direction. However, the ice rule can still be satis- 
fied in low magnetic field. As a result, a reduced degener- 
acy is expected for the kagome plane and the dimension- 
ality of the spin system is reduced. Therefore, the mag- 
netic monopoles can diffuse only in the kagome planed 
as shown in Fig. 1(b). However, it should be noted that 
the experimental investigations on probing the monopole 
excitations have not yet arrived perfect consistency 20 ^ 1 
Studying spin-ice materials by using more different tech- 
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niques is probably an effective way. 

Low-temperature heat transport is a powerful tool to 
probe the properties of elementary excitations 2 ^— and 
the magnetic-field-induced magnetic transitions 3 ^— . In 
principle, the magnetic monopoles, as the elementary ex- 
citations in the spin-ice compounds, can also contribute 
to the heat transport by acting as either heat carriers 
or phonon scatterers. In this work, we study the low-T 
thermal conductivity (re) of DTO single crystals down to 
0.3 K with heat current parallel and perpendicular to the 
(111) plane and in the magnetic field along the [111] axis. 
It is found that the phonon heat transport shows step- 
like decreases at the magnetic phase transitions, which 
is nearly isotropic on the direction of heat current and 
is mainly due to the magnetic scattering on phonons. 
Besides, a remarkable hysteresis of k(H) isotherm ap- 
pears at very low temperatures. The peculiarity is that 
this irreversibility has no direct correspondence with the 
well-known magnetization irreversibility. The possible 
mechanisms are discussed. 



II. EXPERIMENTS 

High-quality DTO single crystals were grown using the 
floating-zone technique. The thermal conductivity at low 
temperatures down to 0.3 K was measured using a con- 
ventional steady-state techniquei 33 ' 34 ' 36 i 37 In this work, 
we measured thermal conductivities along and perpen- 
dicular to the (111) plane, named as Ky and k±, by using 
two samples with sizes of 5.3 x 0.74 x 0.16 mm 3 and 
3.1 x 0.71 x 0.15 mm 3 , respectively. Note that the heat 
currents have to be along the longest dimension, while the 
magnetic field is always perpendicular to the (111) plane 
and is along the shortest and the longest dimension for 
Ky and kj_ samples, respectively. So the demagnetization 
effect is unavoidably different for these two samples. The 
typical misalignment of the magnetic field is less than 2°. 

It is necessary to point out that in this work the 
magnetic-field dependencies of re were measured in a 
static field mode. That is, the measurements were done 
by the following steps: (i) change field slowly to a partic- 
ular value and keep it stable; (2) after the sample tem- 
perature is stabilized, apply a heat power at the free end 
of sample; (iii) wait long enough time until the temper- 
ature gradient on the sample (judged from the time de- 
pendencies of two Ru02 thermometers on the samples) 
is completely stabilized; (iv) record the temperatures of 
two thermometers and obtain the temperature gradient. 
The purpose of this operation is to avoid any misleading 
data. It is known that DTO has very slow spin dynamics 
at low temperatures ^22 which may cause some relaxation 
phenomenon in the temperature gradient measurements, 
if the magnetic excitations are involved in the heat trans- 
port behaviors. Actually, in such thermal conduc- 
tivity measurement in the sweeping-field mode could be 
incorrect because it is impossible to get a steady heat 
flow. 
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FIG. 2: (Color online) Temperature dependencies of the ther- 
mal conductivities of Dy2Ti2 07 in zero field for heat current 
along the (111) plane or perpendicular to it. The inset shows 
the temperature dependencies of the phonon mean free path I 
divided by the averaged sample width W for two sets of data. 



III. RESULTS AND DISCUSSION 

A. Heat carriers 

Figure 2 shows the temperature dependencies of re|| and 
kj_ in zero field. Both kii and k± exhibit clear phonon 
peaks at ~ 15 K and a rough T 2 7 dependence at sub- 
Kelvin temperatures, which is however weaker than the 
standard T 3 behavior of phonon thermal conductivity at 
the boundary scattering limit. 22 Similar result has been 
obtained for heat current along the [110] direction, and 
the magnetic scattering on phonons was discussed to be 
important at temperatures below 10 K.— It is possible to 
estimate the mean free path of phonons at low temper- 
atures and to judge whether the phonons are free from 
microscopic scattering at subKelvin temperatures. The 
phononic thermal conductivity can be expressed by the 
kinetic formula n p h — hCv p l^ where C = /3T 3 is the 
phonon specific heat at low temperatures, v p is the av- 
erage velocity and I is the mean free path of phonons. 
Using the f3 value obtained from specific-heat measure- 
ments (not shown here) , the phonon velocity can be cal- 
culated and then the mean free path is obtained from the 
re i 36 i 41 The inset to Fig. 2 shows the ratios l/W for the 
two samples, where W is the averaged sample widthi 36 i 41 
It is clear that both ratios increase quickly at very low 
temperatures and are expected to approach 1 below 0.3 
K. This estimation tells us that it is not likely that the 
magnetic excitations in DTO (for examples, magnetic 
monopoles) can make a sizable contribution to transport- 
ing heat, since the experimental re is even smaller than 
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FIG. 3: (Color online) Magnetic-field dependencies of the k\\ 
and k± for field along the [111] direction. The data are taken 
with decreasing field from high-field limit (14 T). Panels (a) 
and (c) show the raw data, while (b) and (d) show the data 
with magnetic fields re-scaled by 1.5 and 0.95 T, respectively, 
considering the demagnetization effect. Note that the data 
after this re-scaling are essentially isotropic. 



the phonon term in a boundary scattering limit. It should 
be pointed out that the magnetic-monopole excitations 
were believed to be negligible at very low temperatures, 
because the energy barrier for nipping a spin or creating 
a pair of magnetic monopoles is about 4.4 K.— 



B. k(H) and field-induced magnetic transitions 

Figures 3(a) and 3(c) show the k(H) isotherms for Kn 
and respectively, in magnetic field perpendicular to 
the (111) plane. Note that these data are taken with 
sweeping field down from 14 T. One can easily notice that 
at very low temperatures K\\(H) and n±(H) display two 
anomalies at low magnetic fields. For example, at 0.36 K 
the k±(H) curve shows two clear step- like transitions at 
about 0.25 and 0.9 T, which are in good correspondence 
with the critical fields of the subsequential transitions 
from the low-field spin-ice phase to the kagome-ice phase 
then to the saturated state ("3-in, 1-out" or "3-out, 1- 
in" ) iiSri£ In present work, we denned these two transition 
fields as H c \ and H C 2 . The sharpness of the transitions is 
quickly smeared out upon increasing temperature. Sim- 
ilar results are observed in the low-T Ku(H) isotherms. 
It is clearly seen that the main difference between the km 
and k± is that the critical fields of the step-like transi- 
tions are somewhat different. However, if the horizontal 
axis of the n\\ (H) and n±(H) curves are re-scaled by some 
particular values, as shown in Figs. 3(b) and 3(d), one 
can obtain the essentially same behavior of km and k±. 



This means that the difference between the experimental 
data of Ku and k± is simply due to the demagnetization 
effect^ Therefore, the heat transport behaviors shown 
in Fig. 3 actually indicate a nearly isotropic behavior 
and a close relationship between the heat transport and 
the magnetic phase transitions. 

The drastic change of heat transport at the magnetic 
phase transition has been widely observed in many kinds 
of magnetic materials and is closely associated with the 
evolution of magnetic excitations.— ~— However, the par- 
ticular features of n{H) at the magnetic transitions can 
be significantly different from each other, which depends 
on both the natures of magnetic transitions and the roles 
of magnetic excitations in the heat transport (as heat 
carriers or phonon scatterers). In principle, the mag- 
netic monopoles as a kind of excitations in the spin- 
ice and kagome-ice states can either transport heat or 
scatter phonons. However, it is easy to know that the 
magnetic monopoles are not likely to act as heat carri- 
ers for two reasons. First, as mentioned above, in zero 
field the monopole excitations are actually very difficult 
to emerge at subKelvin temperatures. Second, a low 
field along the [111] direction can significantly enhance 
the excitations of magnetic monopoles, whereas the ther- 
mal conductivity does not show any increase. Then, 
can the step-like decreases of k at two transitions be at- 
tributed to phonon scattering by the sudden increases of 
magnetic monopoles? There is no doubt that the mag- 
netic monopoles can effectively scatter phonons by ab- 
sorbing phonon energy and flipping spins, and suppress 
the phonon heat conductivity. However, above H C 2 (~ 
1-1.5 T), the magnetic monopoles are so populated that 
they fully occupy the kagome network with the monopole 
pairs sitting on the nearest sites, in which case their abil- 
ity of scattering phonons must be negligible. Therefore, 
even if the step-like decrease of k(H) at iJ cl can be ex- 
plained by the enhanced phonon scattering by magnetic 
monopoles, the one at H c2 cannot be related to the mag- 
netic monopoles. 

It is known that, due to the strong Ising anisotropy of 
Dy 3+ spins, the magnetic monopoles connected by the 
Dirac strings are the elementary magnetic excitations of 
the spin ice state. The n(H) behaviors shown in Fig. 
3, however, indicate that there may be other magnetic 
excitations playing an important role in the step-like de- 
creases of k at two magnetic transitions. In this regard, 
there has been no obvious evidence that DTO could have 
other magnetic excitations than the magnetic monopoles, 
like the spin fluctuations caused by non-Ising terms in 
the Hamiltonian i 21 i 43 Nevertheless, some drastic field- 
induced changes related to crystal lattice and phonons 
have been indeed found. In a recent work, the ultra- 
sound measurements also indicated sharp anomalies of 
the sound velocity and sound attenuation at H C 2, which 
may share a common origin with the k(H) transitions.— 
It is also notable that a local maximum of entropy at 
H c 2 may also indicate a strong spin fluctuations that can 
strongly scatter phonons 4^ 
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FIG. 4: (Color online) Ku(H) and k±(H) loops with mag- 
netic field sweeping up and down (perpendicular to the (111) 
plane). The data shown with solid symbols are measured in 
ascending field after the sample is cooled in zero field, while 
the open symbols show the data with descending field, as in- 
dicated by arrows. The demagnetization effects, which are 
the same as those in Fig. 3, are not taken into account in 
these plots. 

Another important impact on the held dependence of 
k that should be taken into account is the paramag- 
netic scattering effect on phonons j 22 i 45 ' 46 It is known that 
the Dy 3+ ions have a degenerate doublet of the lowest 
crystal-held level, which was evidenced as a simple two- 
level Schottky anomaly of the low-T specific heat in mag- 
netic fields! 47 ' 48 Therefore, a Zeeman-effect splitting of 
the doublet can produce resonant scattering on phonons 
and give a low-field suppression of thermal conductivity, 
as some other magnetic materials have showni 45 i 46 This 
can also explain the continuous decrease of k at relatively 
high temperatures (above ~ 1 K), where two magnetic 
transitions are gone. 



C. Irreversibility of k(H) curves 

An unusual phenomenon in DTO is that, as shown in 
Fig. 4, the K\\(H) and k±(H) isotherms show an irre- 
versibility at very low temperatures. For example, after 
cooling the sample to 0.36 K in zero field, both k» and 
k±_ are firstly measured with increasing held, for which 
they show a valley-and-peak feature. Then, upon de- 
creasing field from above 2 T, the magnitudes of k\\ and 
k± deviate from those of the field-increasing curves be- 
low ~ 1.5 and 1 T, respectively. Actually, the differences 
between field-up and field-down curves are so large that 
the valley-and-peak feature of k{H) does not show up at 
all for decreasing field. Upon increasing temperature, the 



irreversibility is quickly diminished and disappears above 
0.7 K. The peculiarity of this transport irreversibility is 
that it has no direct correspondence to that of the low- 
T magnetization curves M(H)^ which feature a pro- 
nounced hysteresis below 0.5 T and a small one around 
0.9 T that were ascribed to the slow dynamics of the 
spin ice and the first-order phase transition, respectively. 
Furthermore, some higher-field hysteresis of M(H) was 
observable only below 0.36 KJ£ 

Although the irreversibility of thermal conductivity 
can sometimes be observed at the field-induced first-order 
phase transitions in many materials, they are usually very 
small and presented in a narrow vicinity of the critical 
fields^ One exception for showing large hysteretic heat 
transport is the high-T c superconductors in the mixed 
state, 49 for which the hysteretic k(H) behaviors are ac- 
companied with the magnetization irreversibility caused 
by the vortex pinning. In other words, the difference 
of k(H) between the field-up curve and the field-down 
one is directly related to the magnetism, which is reason- 
able because the density and distribution of the vortices 
determine the strength of vortex-electron scattering and 
the electron heat transport^ In contrast, the irreversible 
k(H) behaviors of DTO can exist even in the field re- 
gion where the magnetization shows a nearly reversible 
plateau*^ In addition, the ultrasound properties also dis- 
played a broad hysteresis, but only at temperatures as 
low as 0.29 K.— Since the heat transport measurement 
is performed in a static-held mode, in contrast to the re- 
cent interesting nonstationary magnetization and ultra- 
sound measurements, it may reveal different physics 
from other techniques. 

It is necessary to point out that the step-like decreases 
in the low-T k(H) isotherms in some sense are not ex- 
traordinary, since similar anomalies have been observed 
in other magnetic material s 30 ' 37 and are not difficult to 
understand. Although the hysteresis of n(H) looks less 
remarkable than the step-like features in the magnitude, 
it is actually a more peculiar phenomenon related to the 
magnetism of the spin-ice materials. 



D. Possible origins for the irreversibility of k(H) 

One possible origin of the irreversible k(H) behaviors 
is related to field misalignment from the [111] direction, 
which is known to have important impact on the field- 
induced transitions. 51 When the field is slightly tilted 
from the [111] direction due to a misalignment, its com- 
ponent parallel to the (111) plane can break the three- 
fold degeneracy of the modified ice rule. Accordingly, 
the kagome spins tend to align at a temperature that 
depends on this parallel component. This is known as 
a Kasteleyn transition^ 1 - An early entropy measurement 
indicated that a field tilting larger than 2° can signifi- 
cantly smear out the local entropy peak at the H C 2^ 
If such a misalignment of spins could result in a hys- 
teresis in the entropy of the kagome-ice state when the 



field is swept up and down, the phonon heat transport 
would show a hysteresis considering the magnetic scat- 
tering on phonons. Furthermore, the partial ordering 
caused by the parallel component of field is likely ac- 
companied with hysteresis because of the slow spin dy- 
namics. Since this kind of hysteresis may appear only in 
the magnetization component perpendicular to the [111] 
direction, it cannot directly be observed in the magne- 
tization along the [111] direction. Therefore, it would 
not be surprising that the k(H) irreversibility has no di- 
rect correspondence to the well-known M{H) curves. For 
these reasons, the field misalignment should be carefully 
taken into account. However, there are two factors to 
be considered if one tries to ascribe the irreversibility of 
k(H) to the field misalignment. First, since both the 
low-T K\\(H) and k±(H) isotherms in the present work 
show rather sharp decreases at H C 2, the misalignment in 
our measurements might be negligibly small. Second, it 
has been not yet proved that such kind of field misalign- 
ment can really lead to some irreversible behavior in the 
kagome-ice state. 

We propose another possibility with a "pinning" ef- 
fect of the magnetic monopoles by the static disorders, 
in analogy to the weak-disorder induced pinning effect 
on vortices in the superconductors. Considering a sim- 
ple case that the defect is on the nonmagnetic Ti 4+ site, 
the interactions of the Dy 3+ spins in that tetrahedron 
would not be the same as others and the flipping of these 
spins could be more difficult. In other words, a mag- 
netic monopole located at this defect site cannot move 
away freely. Thus, crystal defects can work as pinning 
centers of the magnetic monopoles. Figure 5 illustrates 
the movement process of monopoles considering this kind 
of pinning effect. In zero field, there are few magnetic 
monopoles induced by thermal fluctuations (at very low 
temperatures), as shown in Fig. 5(a). When the field 
is increasing, some spins are flipped and the magnetic 
monopoles are created in pairs (see Fig. 5(b)) and diffuse 
freely in the crystal (see Fig. 5(c)). With increasing the 
field further, the density of magnetic monopoles increases 
and they could contribute to scattering phonons when 
they are moving in the (111) plane; this results in a quick 
suppression of k, as shown in Figs. 4(a) and 4(b). Note 
that in this process, once some magnetic monopoles meet 
with the defects, they can be pinned. Near the saturation 
field, there are so many monopoles that they gradually 
lose their mobility, and the scattering effect on phonons 
becomes weakened and k can be somewhat recovered. 
At the saturated state (the spins are fully polarized at 
H C 2), the magnetic monopoles are well populated so that 
they fully occupy the kagome network with the monopole 
pairs sitting on the nearest sites in Fig. 5(d). At this mo- 
ment, they completely lose the mobility and the k reaches 
the local maximum (the decrease of k at H > H C 2 is 
due to other reasons, as was discussed in section III. B). 
When the field is decreasing across the saturation field, 
the density of magnetic monopoles is quickly reduced by 
annihilation in pairs. Most of remained ones may have 
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FIG. 5: (Color online) Structure evolution of the spin projec- 
tion on the (111) plane. The magnetic field is perpendicular 
to the plane. H c i and H C 2 are defined in the text. Panels 
(a)-(d) show the field-up process while (e)-(f) show the field- 
down process, (a) In zero field and at very low temperatures, 
the thermally excited magnetic monopoles are very few. The 
black circles stand for crystal defects, (b) With increasing 
field, some magnetic monopoles are created in pairs (the red 
and blue dots), (c) Magnetic monopoles can move separately 
by spin flipping. The flipped spins are shown by the open and 
yellow arrows. The dashed red and blue arrows, accompanied 
with the dashed open circles, show the routines of positive and 
negative monopoles. Once some magnetic monopoles meet 
with the defects, they can be pinned, (d) Above H C 2, the 
magnetic monopoles are so well-populated that they cannot 
move any more, (e) With decreasing field, the number of mag- 
netic monopoles is decreased, (f) Those magnetic monopoles 
sitting on the defect sites are pinned, while the others still 
can move in the plane. 



some chance to start moving, but those pinned by the 
defects cannot be mobile, as shown in Figs. 5(e)-(f). As 
a result, the proportion of the mobile monopoles in the 
field-decreasing process is smaller than that in the field- 
increasing case. Consequently, the phonon scattering is 
weaker in the former case and k shows an irreversible 
behavior. It could be expected that the irreversibil- 
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ity of k must be quickly weakened with increasing tem- 
perature, because the stronger thermal fluctuations will 
cause the magnetic monopoles to be depinned. One may 
think that, based on the above discussion, the peculiar 
heat transport properties can be explained by a simple 
phonon-assistant spin flip and the magnetic monopoles 
might be irrelevant. However, it is notable that the irre- 
versibility is larger in ku (H) than that in n±(H). This in- 
dicates that the magnetic monopoles can more effectively 
scatter the in-plane phonons as a kind of quasi-particles 
since they can move in the kagome planed The validity 
of this picture needs further experimental and theoretical 
investigations. 

IV. SUMMARY 

The very-low-temperature thermal conductivity of 
Dy2Ti2 07 displays an irreversibility with magnetic field 
applied in the [111] direction, which seems to have no 
direct correspondence with the magnetization hysteresis. 



We discussed possible origins of this irreversibility, in- 
cluding a field misalignment effect and a pinning effect 
of magnetic monopoles by the weak disorders. 
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Note added. —We note a recent study on the heat trans- 
port of Dy2Ti207, 52 in which the measurements were 
done in magnetic field along the [001] direction. An irre- 
versible behavior of n(H) was also presented in low fields 
(< 0.3 T), which is directly related to the irreversible 
magnetization. This is very different from what we ob- 
served in a much higher and broader field range. 
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